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AERODYNAMIC DERIVATIVES FOR AN OBLIQUE WING AIRCRAFT ESTIMATED 
FROM FLIGHT DATA BY USING A MAXIMUM LIKELIHOOD TECHNIQUE 


Richard E . Maine 
Dryden Flight Research Center 


INTRODUCTION 


The oblique wing configuration is of interest because of its potential for tran- 
sonic drag reduction and corresponding fuel savings. However, there are several 
unanswered questions about the aerodynamic characteristics of oblique wing aircraft 
and the control laws necessary to achieve acceptable handling qualities. Before 
the oblique wing concept can be developed to its full potential it will be necessary to 
develop techniques to verify the aerodynamic characteristics of unsymmetrical 
aircraft in flight. 

The theory of estimation of stability and control derivatives (refs. 1 to 5) 
can be readily applied to oblique wing aircraft by using a five-degree-of-freedom 
model , but computational complexity and difficulties in maneuver selection make 
this approach unattractive. This report describes a technique previously used for 
symmetrical aircraft (ref. 4) that simplifies the estimation by separating the analysis 
of the longitudinal and lateral-directional motions without neglecting the cross- 
coupling terms . The technique was briefly described in reference 6 . 

A flight test program was conducted with a small, remotely piloted, oblique 
wing aircraft . Maneuvers were performed at wing skew angles of 0° to 45° . The 
resulting data were digitally filtered and used to estimate a complete set of stability 
and control derivatives. The derivative estimates obtained are presented and 
compared to predictions . 


SYMBOLS AND ABBREVIATIONS 


All data are referenced to fuselage body axes according to right-handed sign 
conventions . 

a normal acceleration, q 

n ^ 


longitudinal acceleration, g 


lateral acceleration , g 




g( ) 


reference span , m 
dummy variables (p . 6) 

nondimensional lift coefficient 

nondimensional rolling-moment coefficient 

nondimensional pitching-moment coefficient 

nondimensional normal-force coefficient 

nondimensional yawing-moment coefficient 

nondimensional side-force coefficient 

reference chord, m 
general function 

spectral density of measurement noise 

2 

acceleration due to gravity, m/sec 
general function 

2 

moment of inertia about roll axis , kg-m 

2 

propeller moment of inertia about roll axis , kg-m 

2 

cross product of inertia between roll and pitch axes , kg-m 

2 

cross product of inertia between roll and yaw axes, kg-m 

2 

moment of inertia about pitch axis , kg-m 

2 

moment of inertia about yaw axis , kg-m 
cost functional 
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M 

PCM 

PSD 


P 

q 
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T 
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u 

V 

X 

X , X , X , Xa 
a a a B 

n AJ r 


Z 

z 



dummy variables (p • 6) 
angle of attack upwash factor 
mass, kg 

pulse code modulation 
power spectral density 
roll rate, deg/ sec or rad/ sec 
pitch rate, deg/ sec or rad/ sec 

2 

dynamic pressure, N/m 
yaw rate, deg/ sec or rad/ sec 

reference area, m^ 
maneuver duration, sec 
time , sec 
input vector 
velocity, m/sec 
state vector 

longitudinal offsets of instruments from center of gravity , m 

measured observation vector 

discrete time transfer function variable 

vertical offsets of instruments from center of gravity, m 

computed observation vector 


a 

P 

A 


angle of attack , deg 
angle of sideslip , deg 
increment 


H-1003 
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6 
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Superscripts: 

m 

* 


Subscripts: 

c 

P» P> 

»a' «e’ »r 

0 


aileron deflection (left minus right) , deg 

elevator deflection , deg 

rudder deflection, deg 

measurement noise vector 

pitch attitude, deg 

vector of unknown coefficients 

3 

air density , kg/m 
bank angle , deg 
propeller speed, rad/sec 


measured 
matrix adjoint 

derivative with respect to time 
computed observation 

derivative with respect to indicated quantity 
bias 


VEHICLE AND INSTRUMENTATION 


To investigate the feasibility of flying an oblique wing aircraft , the NASA Dryden 
Flight Research Center and the NASA Ames Research Center have conducted flight 
tests with a small oblique wing aircraft, which is referred to hereinafter as the , 
oblique wing (fig. 1 , ref. 7) . The aircraft's wing was capable of being skewed 
up to 45® left wing forward. The aircraft was flown in two configurations, the 
short and long tail versions, with lengths of approximately 6.6 and 7.6 meters, 
respectively (fig. 2) . The geometric characteristics of the vehicle are given in 
table 1. Power was provided by a 90-horsepower, four -cylinder , air-cooled, 
two-stroke, reciprocating engine that drove a ducted midfuselage three-bladed 
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propeller with a cruise speed of approximately 3600 rpm . The propeller rotated 
clockwise viewed from the front of the aircraft . A fixed tricycle landing gear 
allowed horizontal takeoffs and landings on a dry lakebed . 

The aircraft was equipped with an all-movable tail and conventional rudder 
and aileron surfaces that were controlled remotely by a pilot in a ground cockpit. 
Telemetry links were used to send pilot control inputs to the vehicle and to return 
aircraft response data to the ground for the pilot display. In addition, a television 
camera was mounted in the aircraft's nose, and the picture was displayed in 
the ground cockpit. 

The augmentation system consisted of pitch and roll attitude feedback. Three 
attempts to utilize the pitch augmentation resulted in dynamic instabilities, so it was 
not used for the remainder of the program . The roll augmentation was on during 
several of the test maneuvers , but it had no noticeable effect on the estimates . 

The vehicle's weight and center of gravity were measured on balance scales 
for each wing skew angle . The moments of inertia for the long tail configuration 
were estimated from swing tests (ref. 8) of the assembled vehicle after the last 
flight. The moments of inertia for the short tail configuration were estimated 
by subtracting the predicted differences from the long tail values. The mass 
data for the short tail configuration are given in table 2 . For this configuration , 
the fuel tank was near enough to the center of gravity and the fuel weight low 
enough so that the moments of inertia did not change significantly during the data- 
gathering portion of the flight . For the long tail configuration , a second fuel 
tank was added near the front of the vehicle. Consequently, the changes in 
the moments of inertia during the long tail flights were significant . Table 3 
gives the long tail mass data as a function of fuel loading. 

All the standard stability and control variables were Instrumented , including 
static and total pressures , control surface positions , vehicle attitudes and angular 
rates , linear accelerations , and angles of attack and sideslip . The angle of attack 
and sideslip measurements were taken from standard boom-mounted metal vanes, 
and control surface positions were measured by control position transducers on the 
control surfaces . Table 4 shows the vane and accelerometer locations . The loca- 
tions are given relative to the reference center of gravity, which was 0.01 meter 
above the centerline at 29.3 percent of the reference chord. 

The data were sampled at 200 samples per second by an 8-bit pulse code 
modulation (PCM) system and telemetered to the ground for recording. The 
resolutions of the PCM signals are listed in table 5 . 

The angle of attack measurement was unusable for approximately one-third 
of the second flight. Maneuvers obtained during this period were analyzed with- 
out using the angle of attack measurement. No other instrumentation problems 
were observed . 
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DATA PROCESSING 


II ■Hill llllllllllllil 


A digital filter was used after the flight to remove the effects of several 
structural modes that were excited by engine vibration . The filter consisted 
of both a notch at 17.7 hertz and a third-order Butterworth lowpass with a break 
frequency of 20 hertz . 

The z-plane transfer function of the low pass was 

K, (z + 1)^ 

r„ -0.866b^ -1.732b^'l/ -b \ 

z^- 2e ^ cos (0.5bj^)z + e K z - e ) 

where 


b 


L 



and 


i- 


-0.866b, 


2e 


^L = 


The notch transfer function was K 


cos O.Sbj^ + e 






2 cos 


where 


2 

[z^ -2ze 


b =17 7^2L 
°N ^‘-'200 


-0.707b 


+ e } 


and 


K 


N 


4 - 


-0.707b 


2e 


N 


-1.414b 


+ e 


N 


f 2 - 2 cos 


V 




Figure 3 shows power spectral density plots of three representative unfiltered 
signals . The and p plots exhibit several structural modes between 40 hertz and 

80 hertz and one broad mode centered at 17.7 hertz. The strong peak at 60 hertz 
is from engine vibration . A boom bending mode at about 4 . 5 hertz is evident in 
the a signal . The boom mode was not filtered out because its frequency was too 
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low to filter without affecting the aerodynamic modes . The magnitude of the boom 
oscillation was not large enough to cause problems in the analysis . The power 
spectral densities of the filtered data are shown in figure 4. The high frequency 
modes are significantly attenuated . The boom bending mode in the a signal is essen- 
tially unaffected , as expected . Time histories of the longitudinal signals are shown 
before and after filtering in figures 5 and 6, respectively. Marked improvement in 
and q is apparent in the filtered data with no noticeable phase shift or attenuation 

in the aerodynamic responses. The boom structural mode can be observed near the 
end of the filtered a time history . 

If the unfiltered data are thinned below 200 samples per second, aliasing of 
the 40-hertz to 80-hertz structural responses becomes a problem. In the thinned 
data, these responses are folded down to lower frequencies and cannot be separated 
from the rigid body response . The folding is illustrated by figure 7 , which is the 
power spectral density of the unfiltered thinned to 25 samples per second. The 

structural modes from figure 3 have folded down and spread over the spectrum 
between 0 hertz and 12 .5 hertz, with the result that the spectrum is nearly white. 

It would be almost impossible to recover a good signal from these data; therefore, 
200-sample-per-second data were necessary for filtering out the structural response. 
After filtering, the data were thinned to 25 samples per second for analysis. 


FLIGHT TEST PROCEDURE 


The flight program was limited to three flights of approximately 1 hour each. 
Longitudinal stability problems required a configuration change after the first flight 
(the tail was moved back 1 m) , leaving only two flights in the final configuration. 
Since the flights were separated by only a few weeks , only a preliminary data analy- 
sis could be performed between flights. Thus no time was available to experiment 
with different maneuvers for the identification of the stability and control derivatives. 

The pilot performed 87 maneuvers for stability and control data during the three 
flights. Elevator doublets were used for longitudinal maneuvers; each lateral- 
directional maneuver included a rudder doublet and an aileron doublet . The 
instability in the pitch augmentation system caused three unplanned elevator oscilla- 
tions in the first flight which were also used as test maneuvers. Radio interference 
problems resulted in several periods of apparently random control surface motion , 
seven of which were selected as test maneuvers for analysis. Analysis was 
attempted on these 97 maneuvers, some of which involved inputs to both the longitu- 
dinal and lateral-directional control surfaces. 

The phugoid mode was relatively difficult to control because of the small size 
and low speed of the vehicle; thus, many of the maneuvers were obtained when flight 
conditions were not well stabilized. The phugoid problem was worse when the quality 
of the television picture was poor , since the pilot depended heavily on television 
for attitude reference. The radio interference problems mentioned above also 
increased the pilot's workload considerably. 
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PREDICTED AERODYNAMIC CHARACTERISTICS 


Predictions of the static derivatives were based on data obtained from full-scale 
tests conducted in the Ames 40- by 80-Foot Wind Tunnel. The actual flight vehicle 
in the short tail configuration was placed in the wind tunnel to obtain the data. 

Rotary damping derivatives were predicted theoretically by using lifting line theory 
(ref. 9) . Predictions for the long tail configuration were obtained empirically by 
multiplying the short tail predictions of C , C , C , and C by 1.33 (the 

a ^6 6 

^ e r 

2 

ratio of the tail volumes); for C and C a ratio of 1.77 (1.33 ) was used. 

q r 

The wind tunnel data presented herein are corrected for power effects and are 
referenced to a center of gravity of 29.3 percent of the reference chord. 

The vehicle was tested in the wind tunnel in both the short and the long tail 
configurations. There was a significant difference between the wind tunnel and 
flight test versions of the aircraft. The difference concerned the tail boom and the 
truss arrangement that attached the empennage to the vehicle. For the wind tunnel 
tests , the tail boom was composed of several interchangeable parts to facilitate the 
testing of various tail lengths. This caused considerably more tail boom flexibility 
than in the flight vehicle, in which the tail boom was a solid tube and quite rigid. 

The flexibility of the wind tunnel vehicle in the long tail configuration was great 
enough so that the long tail data were judged unusable. Therefore, the simple 
empirical factors mentioned above were used to obtain predicted long tail deri- 
vatives from the short tail predictions. The flexibility of the wind tunnel vehicle 
in the short tail configuration was not as extreme , but it still opens the validity 
of the predictions to question . 

Analysis of the static data was further complicated by the fairly large scatter 
in some of the data. The scatter was due in part to the aircraft's size, which was 
small for the wind tunnel's data acquisition system. For example, over the range 
of angle of attack for which flight data are presented (3° to 9°) , the aircraft's 
pitching-moment coefficient increment is approximately -0.0400 on the basis of the 
static wind tunnel data . This value is only three times or so as large as the quoted 
resolution of each static data point in the range . The scatter was great enough so 
that predictions could vary by a factor of 2 or 3 for the same flight condition 

and configuration , depending on the fairing used and the breakpoints plotted . The 
wind tunnel rolling-moment derivatives also exhibited a large amount of scatter , 
and in some cases the scatter made the predicted roll-related coefficients questionable. 


METHOD OF ANALYSIS 


The maximum likelihood estimation method was used to analyze the oblique 
wing data. This method is widely accepted as one of the best techniques for 
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estimating stability and control derivatives from flight data . The maximum likelihood 
method for systems with measurement noise but no state noise is briefly reviewed 
below . Further discussion is to be found in reference 2 . 

The system is defined by the following equations: 

x(0 = f[x(t) , u(t) , |] 

2(0 =g[x(t), u(t), y +Ti(0 


where 

X state vector 

u input vector 

^ vector of unknown coefficients 

z measured observation vector 

r| measurement noise vector 

If T) is assumed to be band limited, white, zero mean, Gaussian noise with spectral 
density G , the maximum likelihood estimate is obtained by choosing ^ to minimize the 
cost functional 


J(V 


= /o"[ 


z ^(0 - 2(0 


ij*(GG*) 


L(t) - z(0 \dt 


where z^ is the observation vector obtained by integrating the system equations with 

the assumed value of ^ and no noise. An extra term is often included in maximum 
likelihood estimation routines to include the effects of a priori information , but this 
feature was not used for the analysis described in this paper. 


A modified Newton-Raphson (Newton-Balakrishnan) iterative algorithm is used 
to minimize J (refs . 2 and 3) . 


Five-Degree-of-Freedom Approach 

The five-degree-of- freedom approach to the analysis of oblique wing data is 
the most straightforward approach, and it is suggested by many investigators. 
However, because of the difficulties inherent in the five-degree-of-freedom approach, 
it was not used on the oblique wing data; instead, a simpler approach was sought. 

It is generally accepted that five-degree-of-freedom equations of motion are 
necessary for the analysis of oblique wing data (the sixth degree of freedom , velocity , 
can still be ignored) . Three types of coupling between the longitudinal and lateral- 
directional motions combine to invalidate the standard three-degree-of-freedom 
assumptions . The kinematic coupling terms present but often ignored for symmetrical 
aircraft are larger for the oblique wing, because for the oblique wing, significant 
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motion is present in all axes. Kinematic coupling terms not present for symmetrical 
aircraft arise because of the nonzero I^y of the oblique wing. Finally, the oblique 

wing results in significant aerodynamic coupling between the modes . 


The application of the maximum likelihood estimation method to the five-degree- 
of-freedom equations for the oblique wing is straightforward from a theoretical 
standpoint; however, problems abound in the practical application of the technique. 
There are seven states, nine observations, and three controls. With the derivatives 
we choose to estimate there would be 39 unknowns , so computation time would be 
relatively long. Some numerical problems are also to be expected from handling 
equations with 39 unknowns . Although this number of unknowns is not large 
enough to make the problems extreme , double precision would almost surely be 
required if a 32-bit computer were used; naturally, double precision would not help 
the computation time problem . 


Maneuvers must be carefully selected and performed to excite the modes enough 
to permit an accurate simultaneous estimation of all the coefficients . The selection 
and refinement of the maneuvers would be particularly difficult in a short flight 
program such as that for the oblique wing . 

Overall , it was judged that it would be difficult to obtain estimates of high quality 
with the five-degree-of-freedom approach. The inaccuracy of the linear aerodynamic 
model, the brevity of the flight program, and the difficulties of the five-degree-of- 
freedom approach combined to make such an analysis a major effort. Therefore, 
a simpler technique was sought , and the five-degree-of-freedom approach was not 
used on the oblique wing data . 


Analysis by Separation of Modes 

Many of the problems arising from the five-degree-of-freedom approach could 
be alleviated if the longitudinal and lateral-directional motions could be analyzed 
separately . For conventional aircraft , _the motions are usually separated by simply 
neglecting the cross-coupling terms. For an oblique wing, these terms are too 
large to neglect. 

However , for this aircraft , it is possible to include the effect of the cross- 
coupling terms without using all five degrees of freedom . This is done in the 
longitudinal analysis by eliminating the differential equations for the lateral- 
directional motion and using the measured lateral-directional responses as inputs to 
the longitudinal equations. Similarly, the lateral-directional analysis uses the 
measured longitudinal responses. 

The longitudinal equations of motion for the oblique wing , including propeller 
gyroscopic effects , then become 

6 = -A^c. + q + ft (cos e cos cp"* cos a + sin 0 sin a) - tan P™(p'" cos a + r'" sin a) 
mV L, V ' ' ^ ' 


qly -pi 


XY = 


, m m,j j 

+ r p i.ly-1. 


) + 


[ . m.2 , 

(r ) - (p ) J 


XZ 


+ qr 


. ^ r 

Vx, 


• m m . fn 

0 = q cos <p - sin tp 
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The longitudinal aerodynamic parameters chosen to be estimated for the oblique 
wing are contained in the following expansions of the nondimen sional moments 
and forces . The reasons for this choice are explained in the section entitled 
Results and Discussion , 


Cj =Cj a + Cr 6 + C. + C, P 
e 


a + C C 6+C +C 6+C + C + C p'” 

'"/,2V m- e rrin "Ir a 2V 2V 

a q d 06 p r B 

^ e a ^ ^ 


^( 1 ) 


The longitudinal observation equations are: 


a 


c 


= K a 
a> 





0=0 

c 




It should be noted that although the d equation above contains (lift) derivatives, 

this report presents (normal-force) derivatives because the sixth degree of 

freedom, velocity, is not included. For the low angles of attack investigated, this 
substitution should not cause problems . 

The lateral-directional equations are as follows: 

P = cos sin cp + p sin - r cos 

cp = p + r cos cp tan 0^^ + sin cp tan 0^ 
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where the following coefficient expansions contain the lateral-directional aero- 
dynamic parameters chosen to be estimated. 


Cy - Cy P + Cy 6^ + Cy 8^ + Cy + Cy O 
Bo o 0 a 


, , in 

p p r 6^ 8^ 0 q a 

^ , m 

* %W • • V 


y (2) 
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The lateral-directional observation equations are: 

Pe'P-^P^*- 


Pc'P 


r = r 
c 


CPc = 9 


a = ^ p + — ^ r + a 

^9 y g ^ 9 y 
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In the equations above, the m superscripts denote measured quantities. Both the 
longitudinal and lateral-directional equations assume small angles of sideslip . 

Two assumptions are necessary for this approach. The most important assump- 
tion is that measurements of all of the state variables of the unmodeled modes are 
available and contain relatively low levels of noise contamination. The second 
assumption is that the magnitudes of the cross-coupling terms are small compared 
with the standard terms. The second assumption ensures that the errors in the 
cross-coupling terms caused by using the noise-contaminated measurements do 
not significantly affect the analysis. Obviously, there is a tradeoff between these 
two assumptions: The lower the measurement noise, the larger the allowable 
magnitude of the coupling terms . 

This separation technique has previously been used to account for kinematic 
coupling in the symmetrical aircraft (ref. 4); however, the assumptions inherent 
in the technique have not been carefully documented . Examination of the assump- 
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tions shows that the application of the separation technique to the oblique wing 
requires much more caution than the application of the technique to symmetrical 
aircraft because of the increased magnitude of the cross-coupling terms. The 
allowable magnitude of the measurement noise is, therefore, much smaller for the 
oblique wing . 

A slight extension of this separation technique which allows the system to be 
represented in a linearized manner is accomplished by using measured data from 
both the longitudinal and lateral-directional modes in the nonlinear terms. The 
data presented in this report were obtained from equations linearized in this 
manner by using the following specific procedure. The nonlinear terms were simply 
computed by using measured data, with two exceptions. The cos 0 in the a 

equation was replaced by cos 0^ - (0 - 0^) sin 0^ and the sin cp in the p equation 

was replaced by sin + (q) - cp^) cos cp^. This significantly reduces the computer 
time required. It also allows a standard linear program to be used for the derivative 
estimation . 

The computer program used , MMLE 3 , is an outgrowth of the MMLE program 
(ref. 3) . The MMLE 3 program was not written with the oblique wing application 
in mind , but rather to provide a standard identification program with a general 
set of linear equations . 

It should be possible to incorporate the separation technique proposed into 
any reasonably general maximum likelihood estimation program . 


RESULTS AND DISCUSSION 


Coupling between the longitudinal and lateral-directional modes was significant 
on the oblique wing , as expected . Figure 8 shows a typical time history of the 
aircraft's response to an elevator input; motion in all of the modes is considerable 
because of the coupling. Similarly, figure 9 shows that the aircraft responds in both 
the longitudinal and lateral-directional modes to a maneuver involving rudder and 
aileron doublets . 

The data were analyzed first by ignoring all of the coupling terms. The values 
used for the G matrix are given in table 6 . Although this produced reasonable 
results for most of the maneuvers at the lower wing skew angles, it was clearly 
unacceptable in several cases, particularly for those at 30° and 45° of wing skew. 
One of the unacceptable fits is shown in figure 10. 

Then the separation technique was used to include cross-coupling effects in 
the analysis. The G matrix values given in table 6 were also used for this analysis. 
Since aerodynamic, as well as kinematic, coupling was significant, it was first 
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necessary to choose the set of derivatives to include in the model . Considerable 
effort was expended in making this choice. For the final analysis, the derivatives 
in equations (1) and (2) were chosen. Analysis was Started with kinematic 
coupling only , and derivatives were added to the model one at a time to evaluate 
their effects. The choice of derivatives was based on physical reasoning, the size 
of the predicted effects , and the improvement observed in the analysis when each 
derivative was added to the model. 

With these aerodynamic terms included , all 97 maneuvers obtained were 
successfully analyzed. This 100-percent utilization was quite unexpected for this 
unconventional aircraft. Figure 11 shows a resulting lateral-directional fit. The 
fit is good compared with figure 10, which is the fit of the same maneuver with the 
cross-coupling terms omitted. Figure 12 shows a longitudinal fit, which is 
excellent . 

The final derivative estimates are presented in figures 13 to 20. The vertical 
bars in these figures are uncertainty levels (refs. 4 and 5) , which indicate the 
approximate accuracy that can be expected of the estimates. The flight estimates 
in these figures are plotted against an angle of attack adjusted for upwash by 
dividing the average measured angle of attack for each maneuver by a factor of 
1.25. This factor was obtained by averaging the upwash factors estimated for all 
of the longitudinal maneuvers . 

The estimate of the upwash factor for each maneuver was obtained by treating 
it as one of the unknown parameters for each longitudinal maneuver. No flight 
estimate of the bias in the angle of attack was available . 

The derivative C is assumed to be zero in the flight estimation model; 

m m 

(X 

thus, the flight- estimated C includes the effect of any actual C . The predicted 

^ m m. 

q a 

values of C + C are, therefore, used to compare with flight-estimated C 
m. PS 

e sc 

^ 4aT^L used for comparison 

with the flight-estimated . For simplicity, the correction to the predicted 

a 

C was made by using a nominal mass of 400 kilograms and an air density of 

3 

1.104 kg/m , which corresponds to a density altitude of 1067 meters. 


Short Tail Derivatives 

Figure 13 shows the estimates of the lateral-directional derivatives for the short 
tail configuration. Flight data in this configuration were obtained at 0® and 15® of 
wing skew angle. The estimates are in reasonably good agreement with predictions, 
with the notable exception of the rudder derivatives. The flight estimates of 
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C are considerably larger in magnitude than the wind tunnel estimates. The 
"6 

r 

flight data also show an increase in rudder effectiveness as angle of attack increases; 
this trend is not present in the wind tunnel data. The same discrepancies can be 
observed in C„ , although the wind tunnel estimates lie within the uncertainty 
6 

r 

envelope. The tail flexibility during the wind tunnel tests may account for some of 
these discrepancies. The flight estimates of are opposite in sign from the wind 

8 

r 

tunnel estimates . It would be natural to expect a positive value of Cp for a 

r 

rudder above the fuselage; this is in agreement with the flight estimates. The 
sign conventions used for the wind tunnel data were checked for consistency with 
those used for this report, but no discrepancy was found. 


Since the short tail flight data were obtained only at wing skew angles of 0® and 

15°, the effects of skew angle were not large. The effects of skew , which can be 

discerned in Cg and Cg , are in excellent agreement with the predicted effects of 

P 8 

^ a 

skew on these derivatives. 


The lateral-directional cross-coupling derivatives due to longitudinal motion 
are presented for the short tail configuration in figure 14. The flight estimates of 
are considerably larger in magnitude than predicted, although Cg is smaller 

"q q 

than predicted . The flight estimates of Cg and C may be partially compensating 

q \ 

for a modeling error . The angle of attack derivatives are in good agreement with 
predictions; the large apparent scatter in Cg is due to its small magnitude for 

a 


0° and 15® of wing skew. 


Figure 15 summarizes the short tail longitudinal derivatives. The estimates 
of Cj^ are slightly larger than predicted. The C estimates are in good agreement 

with the wind tunnel fairing presented here, although, as discussed above, the 
predicted values of C varied widely. 


The flight estimates of C are larger in magnitude than predicted, perhaps 

mg 

e 

because of the tail flexibility during the wind tunnel tests. The larger than predicted 
elevator effectiveness was a major contributor to a longitudinal sensitivity problem 
that almost caused the aircraft to crash during the first flight. The flight estimates 
of C are somewhat larger in magnitude than the predictions. 

q 
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Figure 16 presents the short tail data for the longitudinal cross-coupling 
derivatives due to lateral-directional motion . The agreement with predictions 
is good for all of these derivatives. 

Preliminary derivative estimates were obtained approximately 2 weeks after 
the first flight. These estimates were used to verify the source of the longitudinal 
sensitivity problem and to update the flight simulator . The pilot commented that 
the simulator flew much more like the aircraft after being updated. The consis- 
tent trends in the final short tail estimates, their low uncertainty levels, and the 
good fits all indicate that the estimates obtained are reasonably good . The 
simulator experience provides additional verification of the quality of the estimates . 


Long Tail Derivatives 

Figure 17 shows the long tail estimates for the lateral-directional derivatives. 
Long tail flight data were gathered at wing skew angles up to 45®; thus, the skew 
effects were more easily observed than in the short tail data . The derivative 
C becomes slightly smaller and C„ more negative as the wing skew angle 

P y 

increases . The derivatives and Cg are smaller in magnitude at the larger 

p 

skew angles , as expected . No definite effect of skew angle is apparent in 

P 

or C . The derivative Cg agrees well with the predictions, decreasing at 

a 

the larger skew angles. The C estimates at 0® and 15® of wing skew average 

"6 

a 

about 0.0001 for the long tail flight data, as opposed to 0 for the short tail flight 
data. This difference is almost surely due to a slight error in the estimated 
inertias for one or both configurations . It is difficult to imagine a significant 
effect of tail position on the actual C ; on the other hand , the C estimate 

a a 

from flight data is very sensitive to (for any estimation technique) . The 
flight C from the long tail data shows a peak at 30® of wing skew , in agree- 
a 

ment with predicted trends , although the magnitudes of the flight and predicted 
values are somewhat different. The rudder derivatives show the same disagreement 
with predictions as was noted in the short tail data . In addition , the long tail 
data show that rudder effectiveness decreases as wing skew angle increases. 

Figure 18 gives the long tail lateral-directional cross-coupling derivatives 
due to longitudinal motion. As in the short tail data, the agreement with pre- 
dictions is poor for the q derivatives and good for the a derivatives. 

The estimates of the longitudinal derivatives for the long tail configuration 
are presented in figure 19. Both and become smaller in magnitude as 
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wing skew angle increases. The disagreement of the flight and wind tunnel 
estimates of C is not surprising , considering the resolution problems with the 

wind tunnel test data for the estimation of C .In addition , the factor used to 

get long tail predicted values from the short tail values (1.33) is not very accurate 
when applied to ; it neglects the change in downwash and the wing pitch- 
ing moment. As in the short tail data, and are larger in magnitude 

6 q 

e ^ 

than predicted . The C estimates are more negative at the larger skew angles . 

q 

Figure 20 shows the long tail longitudinal cross-coupling derivatives due to 
lateral-directional motion . The estimates of C , C„ , and C are in good 

nip iVp mg^ 

agreement with predictions. The derivative C is less negative than predicted. 

mp 

The uncertainty levels and scatter of C are fairly large, although at least its 

sign seems well established . This indicates that C is probably the weakest 

of the cross-coupling derivatives estimated . 

The trends of the short tail estimates were substantiated by the long tail data. 
In addition, the long tail data included maneuvers at up to 45° of wing skew angle. 
As a result, several trends were evident in the long tail estimates that were not 
observed in the short tail data . 


CONCLUDING REMARKS 


A technique was described that estimates aircraft aerodynamic derivatives by 
separating the analysis of longitudinal and lateral-directional motion even in the 
presence of the moderate cross coupling characteristic of oblique wing aircraft. The 
technique has several advantages over the usual five-degree-of-freedom approach 
for some situations. Although previously applied to symmetrical aircraft, this 
technique was not expected to be adequate for oblique wing vehicles. 

This separation technique was applied to obtain a complete set of stability and 
control derivative estimates from flight data for a small remotely piloted oblique wing 
aircraft. This application represents the first time that aerodynamic derivatives have 
been estimated from flight data for such aircraft. Aliasing problems necessitated a 
200-sample-per-second rate for the raw data, even though a 25-sample-per-second 
rate for the filtered data was sufficient for the analysis. The complete set of flight- 
determined estimates was presented and compared to predictions. Comparisons with 
predictions were generally reasonable, with some exceptions. The limited number 
of flights available and the lack of usable wind tunnel data for one of the aircraft 
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configurations prevented the comparison from being definitive . However , the 
results demonstrated that the relatively simple approach developed was adequate 
to obtain usable estimates of the aerodynamic derivatives of oblique wing aircraft. 


Dry den Flight Research Center 
National Aeronautics and Space Administration 
Edwards, Calif., December 6, 1977 
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TABLE 1.— GEOMETRIC CHARACTERISTICS OF OBLIQUE WING 


2 

Reference area, m 9.27 

Reference (root) chord, m 1.73 

Reference span, 6.81 

Horizontal stabilizer, m— 

Span 2.68 

Root chord 0.686 

Vertical tail, m— 

Span 0.65 

Root chord 0.94 

Tip chord 0.40 

Rudder span, 0.45 

Aileron span, 1.83 

Propeller diameter, m 1.23 

Reference center of gravity— 

Percent of reference chord 29.3 

Distance behind nose tip , m 2.286 

Vertical center of gravity, 

distance above centerline, m 0.010 (±0.005) 


TABLE 2. -MASS DATA FOR OBLIQUE WING SHORT TAIL CONFIGURATION 

[Full fuel condition] 

(a) Constant data 


Mass , kg 388 

Center of gravity , 

percent of reference chord . 29.3 

Iv » kg-m^ ............... 0.241 

Ap 


(b) Data dependent on wing skew angle 



Wing skew angle, deg 


0 

15 

r 1 2 

kg-m 

314 

302 

r 1 2 

ly, kg-m 

488 

510 

r 1 2 

kg-m 

698 

698 

r 1 2 

^XZ’ 

69 

69 

r 1 2 

I^Y ’ 

0 

~54 
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TABLE 3. -MASS DATA FOR OBLIQUE WING LONG TAIL CONFIGURATION 


(a) Empty weight data 



Wing skew angle, deg 

0 

15 

30 

45 

r 1 2 

I^, kg-m 

314 

302 

262 

211 

r 1 2 

ly, kg-m 

602 

622 

648 

718 

r 1 2 

I^, kg-m 

824 

824 

824 

824 

r 1 2 

Ixz> 

69 

69 

69 

69 

r t 2 

IxY’ 

0 

-54 

-94 

-108 

r 1 2 

I„ , kg-m 

Ap 

0.241 

0.241 

0.241 

0.241 


(b) Fuel loading data 


Time, 

Event 

('q \ 

Total mass, 

Center of gravity, 

and , 

^>Af 

‘XZ’ 

min 


kg 

percent c 

1 2 
kg-m 

, 2 
kg-m 

0 

Start of takeoff roll 

414 

28.0 

54 

3 

36 

Forward fuel expended 

400 

31.5 

12 

3 

78 

All fuel expended 

383 

33,2 

3 

0 


^Assumes a fuel flow rate of 0.39 kg/min. 

^The values in these columns are added to the empty weight data 
in table 3(a) to obtain the data adjusted for fuel loading. 




O TX) R 


TABLE 4. -INSTRUMENT LOCATIONS RELATIVE TO REFERENCE CENTER OF GRAVITY 


Measured 

quantity 

Distance of instrument 

Distance of instrument 

forward of reference 

below reference 

center of gravity, m 

center of gravity, m 

a . . . . 

2.819 


p . . . . 

2.730 

0.305 

“n • • • 

0.584 


a ... 

X 


-0.098 

a ... 

0.651 

-0.098 

y 



TABLE 5. -PCM SIGNAL RESOLUTIONS FOR OBLIQUE WING DATA 


! 

Signal 

Resolution 

a, deg 

0.15 

p. deg 

0.15 

p , deg/ sec 

0.4 

q , deg/ sec 

0.4 

r, deg/ sec 

0.2 

Q , deg 

0.6 

cp , deg 

0.7 

9 

0.02 

^x’ 9 

0.01 

^y> 9 

0.01 

Sq. deg 

0.1 

6g. deg 

0.1 

6^. deg 

0.15 

V , m/ sec 

0.25 

Altitude , m 

10 
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TABLE 6.-G MATRIX VALUES USED IN ANALYZING OBLIQUE WING DATA. 


Signal 

^^^Diagonal element 
of G~^ 

a , per deg 

40 

q, per deg/ sec 

24 

0 , per deg 

34 

a„. per g 

380 

P , per deg 

20 

p, per deg/ sec 

12 

r, per deg/ sec 

15 

9 , per deg 

13 

a , per g 

650 


Off-diagonal elements are all 0. 





Figure 1. Oblique wing . 
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Tail location for 
flights 2 and 3 
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10-71 1 1 1 1 

10‘^ lo'^ 10° 10^ 10^ 

Frequency, Hz 

Figure 3, Power spectral densities of unfiltered 200-sample-per -second 
data. 
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Frequency, Hz 


Figure 4. Power spectral densities of filtered 200-sample-per-second 
data. 








Frequency, Hz 


Figure 7. Power spectral density of unfiltered 25-sample-per-second 
normal acceleration. 
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Figure 11. Comparison of measured and computed 
lateral-directional motions of oblique wing 
aircraft with 45° of wing skew and cross -coupling 
terms included. 
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Figure 12. Comparison of measured and computed 
longitudinal motions of oblique wing aircraft 
with 45° of wing skew and cross -coupling terms 
included. 
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Figure 13. Lateral-directional derivatives of oblique wing aircraft in 
short tail configuration. 
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Figure 13. Continued. 





Sk6W 

Flight Prediction angle, 
““ deg 

0 » e 0 

m B a 15 

1 Uncertainty level 


0.0000 


- 0.0002 


c -0.0004 

r 

per deg .q.qoob 


-0.0008 


- 0.0010 


$ 


[D|Il 


It (|l 




A. 


6 

a, deg 


8 9 


Figure 13. Concluded. 
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Figure 14. Continued. 
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Figure 15. Longitudinal derivatives of oblique wing aircraft in short 
tail configuration. 
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Figure 15 
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Figure 15. Concluded. 
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Figure 17. Lateral-directional derivatives of oblique wing aircraft 
in long tail configuration. 
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Figure 17. Continued. 
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Figure 17. Concluded. 
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Figure 18. Lateral- directional cross-coupling derivatives of oblique wing 
aircraft in long tail configuration. 
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Figure 18. Continued. 
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Figure 19. Longitudinal derivatives of oblique wing aircraft 
in long tail configuration. 
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Figure 20. Longitudinal cross-coupling derivatives of oblique wing 
aircraft in long tail configuration. 
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Figure 20. Concluded. 
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